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We use the non-relativistic random phase approximation with exchange to perform calculations 
of valence shell photoionization of Ne, Ar, Kr and Xe from their respective thresholds to photon 
energy of 200 eV. The energy derivative of the complex phase of the photoionization matrix elements 
is converted to the photoelectron group delay that can be measured in attosecond streaking or two- 
photon sideband interference experiments. Comparison with reported time delay measurements in 
Ne and Ar at a few selected photon energies is made. Systematic mapping of time delay across 
a wide range of photon energies in several atomic targets allows to highlight important aspects of 
fundamental atomic physics that can be probed by attosecond time delay measurements. 
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I. INTRODUCTION 

Time delay in atomic photoionization has become an 
active and rapidly expanding field of research following 
pioneering experiments on attosecond streaking [1] and 
two-photon sideband interference [2]. Both techniques 
use the XUV pump pulse to ionize the target atom and 
the IR probe to obtain the timing information on the pho- 
toemission process. In attosecond streaking, the vary- 
ing time delay between the the pump and probe pulses 
is mapped onto the photoelectron kinetic energy. The 
whole valence shell is projected onto a photoelectron ki- 
netic energy map (the so-called spectrogram) which is 
then modeled, in the strong field or Coulomb- Volkov ap- 
proximations, with the photoionization time delay being 
treated as a fitting parameter. This measurement re- 
vealed a relative time delay of 21 ± 5 as between photoe- 
mission from the 2p and 2s sub-shells in Ne at 106 eV 
photon energy. The positive sign of the relative time de- 
lay indicates that emission of the photoelectron from the 
2p sub-shell is seemingly delayed relative to that from the 
2s sub-shell. 

In the two-photon intcrferomctric technique, the vary- 
ing time delay between the pump and probe pulses is 
mapped onto the two-photon sideband (SB) oscillations. 
The phase of these oscillations depends on the phase dif- 
ference of the two neighboring harmonics and the time 
delay in atomic photoionization process. The atomic 
time delay can be presented as the sum of time delays 
in the XUV photon absorption and subsequent IR pho- 
ton absorption (continuum-continuum or CC transition). 



ta — Tw + tcc 



(1) 



The Tw term represents the Eiscnbud-Wigncr-Smith 
time delay (or Wigner time delay or photoelectron group 
delay, all these terms will be used interchangeably in the 
present context) which is defined as the energy deriva- 
tive of the complex phase of the quantum amplitude of 
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XUV absorption [1, 3]. More details on the Wigner time 
delay theory can be found in the review article [4]. The 
tcc term is modeled using the lowest order perturba- 
tion theory and asymptotic forms of the continuum wave 
functions thus allowing to obtain the former from an ex- 
perimental measurement [5] . By reconstructing the oscil- 
lations of the SB 22 to 26 of the titanium:sapphire laser 
at 800 nm, Kliinder et al [2] reported the relative time 
delay between the photoelectron emission from the 3s 
and 3p sub-shells of Ar in the photon energy range of 34 
to 40 eV. Whether the 3p electron was delayed relative 
to the 3s one or vice versa was found to depend on the 
photon energy. This measurement was repeated later by 
Guenot et al. [6] and the sign of the relative time delay 
was reverted with the 3s photoelectron being delayed rel- 
ative to the 3p one near the top end of the photon energy 
scale. 

This repeated measurement was prompted by obser- 
vation that the photon energy of 40 eV fell very close 
to the Cooper minimum of the 3s shell. Photoionization 
process in this region is driven very strongly by the many- 
electron correlation between the 3s and 3p sub-shells [7] . 
Such a process cannot be theoretically described using an 
independent electron model like the Hartree-Fock (HF) 
theory. So the interpretation of the two-photon interfer- 
ometric measurement [2] based on this theory had to be 
re-evaluated. A more adequate model that accounts for 
inter-shell correlation in noble gas atoms is the random 
phase approximation with exchange (RPAE or, shorter, 
RPA, both acronyms are used here interchangeably) [8]. 
However, even after including the RPA corrections, the 
agreement between theory and experiment did not im- 
prove [6]. 

Theoretical interpretation of the attosecond streaking 
measurement of Schultzc et al [1] is also not straightfor- 
ward. The group delay difference between the 2p and 2s 
sub-shells in Ne calculated in the HF approximation is 
only 6.2 as [3]. With the added RPA correction of 2.2 as, 
it accounts for less that a half of the experimental value 
of 21 ± 5 as. More accurate simulations that accounted 
for both the XUV and IR fields returned somewhat larger 
values of 10.2 ± 1.3 as [9] and ~ 12 as [10]. These values 
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are still far too small to match the experimental result. 

Even though the streaking IR field is relatively weak, 
its interplay with the long-range Coulomb potential of the 
ionic core (the so-called Coulomb-laser coupling - CLC) 
makes an additional contribution to the streaking time 
delay [11]-[14]. Similar to Eq. (1), the streaking time 
delay can be written as 



T~W + 1~CLC 



(2) 



It was suggested in Ref. [11] that tclc should also in- 
clude the effect of the short-range part of the core po- 
tential and hence Eq. (2) should be modified to contain 
twice the Wigner time delay. This would have resolved 
the difference between the theoretical [1, 3] and experi- 
mental [1] time delays. However, subsequent investiga- 
tion on the model two-electron system, that mimicked 
the energy levels of the valence shell of Ne, proved that 
this conjecture is invalid and that Eq. (2) holds. There- 
fore the controversy surrounding the experiment [1] still 
remains unresolved. 

In the present paper, we concentrate on the Wigner 
component Tyy which enters the atomic time delay (1) 
and the streaking time delay (2) measured in the at- 
tosecond interferometric and streaking experiments, re- 
spectively. The corresponding corrections to the Wigner 
time delay, tqc an d T CLC are more or less universal and 
can be readily evaluated [10, 14]. For two-electron atomic 
transitions, like photoionization with excitation and dou- 
ble photoionization, which arc strongly driven by electron 
correlation, the streaking time delay (2) is further mod- 
ified by the CLC effect on the inter-electron interaction 
[15]. These two-electron processes, however, are outside 
the scope of the present study. The target polarization 
by the streaking IR field can also be safely ignored as it 
should be minimal for tightly bound closed shell atoms. 

In the present work, we perform systematic investiga- 
tion of the Wigner time delay in a series of noble gas 
atoms from Ne to Xe across a wide range of photon en- 
ergies. We demonstrate that in heavier noble gases, be- 
yond Ne, the inter-shell correlation, in the form of di- 
rect Coulomb interaction between atomic electrons as- 
signed to different valence sub-shells, has a strong effect 
on photoionization process in general, and the Wigner 
time delay, in particular. To account for this direct inter- 
electron interaction, we employ the RPA method [16]. 
This method can be viewed as an extension of the HF the- 
ory. The latter accounts for the Coulomb inter-electron 
interaction only indirectly by including some part of it 
in the self-consistent one-electron potential. On the con- 
trary, the RPA method accounts for a significant part of 
the direct inter-electron interaction that results in cre- 
ation of pairwise electron-hole excitations. When more 
complex excitations of two-electron-two- whole states are 
important (see e.g. [17]), alternative methods like R- 
matrix [18] can provide more accurate results. 

We validate our computational technique by making 
an extensive comparison between the calculated and ex- 
perimental valence shell photoionization cross-sections. 



Based on this validation, we make specific predictions 
for the Wigner time delay and perform further compar- 
ison with available experimental time delay data. More 
generally, we demonstrate that the Wigner time delay 
contains important phase information that enables at- 
tosecond time delay measurements to reveal various fun- 
damental aspects of atomic physics. 

The paper is organized as follows. In Sec. II we in- 
troduce our computational models for the independent 
electron descriptions and that with account for the inter- 
shell correlations. In Sec. Ill we present our numerical 
results for outer valence ns and np sub-shells in Ne and 
Ar and ns, np, (n — l)d sub-shells in Kr and Xe. We 
conclude in Sec. IV by revealing the systematic trends in 
time delay of noble gases driven by the peculiarities of the 
clastic scattering phases and many-electron correlations. 



II. THEORETICAL MODEL 

1. Independent electron HF model 

We adopt the photoionization formalism as outlined 
in the monograph [16]. We evaluate the one-photon 
dipole matrix element (^1 ^|-S|</>i) of the transition from 
a bound state i to an incoming continuous wave with 
the given photoelectron momentum k. The magnitude 
of the momentum is restricted by the energy conserva- 
tion E = k 2 /2 = uj + Si , where u> is the photon energy. 
The atomic units are used throughout the paper with 
e = m = h = 1 and the atomic unit of time is approxi- 
mately equal to 24.2 as. 

We split the radial and angular dependence in the ini- 
tial state <t>i{r) = Yi imi ('r)R ni i i {r) and use the partial 
wave expansion in the final state 



3/2 



£1/2 



^^^YLiWU^Rmir) ,(3) 



im 



where the radial orbitals are normalized to energy 
(El\\E'l) = S(E — E') and have the asymptotics at in- 
finity 
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We align the quantization axis z with the polarization 
axis of light and write the dipole operator in the length 
gauge as z — ^/47r/3rYio(r) . We perform the spherical 
integration to arrive to the following expression: 
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Here the reduced dipole matrix element, stripped of all 
the angular momentum projections, is defined as 

(El\\r\\n i l i )=U i ( l Q Jo)/ r 2 dr R E i(r)r ^(r) , 

where I = y/2l + 1. The partial photoionization cross 
section for the transition from an occupied state n,Z, to 
the photoelectron continuum state El is calculated as 

v ni ii->Ei{u) = ^n 2 aaluj \(El \\ r Wriik)] 2 . (6) 

Here a is the fine structure constant and ao is the Bohr 
radius. 

The basis of occupied atomic states ||n^) is defined 
by the self-consistent HF method and calculated using 
the computer code [19]. The continuum electron orbitals 
(El\\ are defined within the frozen-core HF approxima- 
tion and evaluated using the computer code [20]. These 
states are found in the combined field of the nucleus and 
the HF potential of the frozen electron core. So the 
photoelectron scattering phase Si(E) delivered by this 
method contains both the long-range Coulomb and the 
short-range Hartree-Fock components. 

We note that the reduced matrix element (5) is real and 
thus the complex phase of the dipole matrix element (4) 
is defined by the scattering phases Si i ±i(E). According 
to the Fano's propensity rule [21], the dipole transition 
with the increased momentum I = Z, + 1 is usually domi- 
nant. In such a situation, the photoemission group delay 
is approximately given by tw = dSi/dE. 

2. Inter-shell correlation 



2U — W. That explains the term exchange in the name 
RPA(E). The direct Coulomb matrix is expressed as 

(nAElWUWE'l'rHk = W 'tj, ^ q ) ( ) 
x R$ Mi (E,E',n i ,n j ) , (8) 

where is a Slater integral [16]. In the exchange ma- 
trix, the electron El and the hole rijlj states are swapped. 

The RPA equations are represented graphically in Fig. 
1. Here the straight line with an arrow to the left or right 
represents electron (continuum) or hole (bound) states, 
respectively. The wavy line exhibits the Coulomb inter- 
action. The dashed line is used to display a photon of the 
frequency ui. The shaded circle is used to represent the 
correlated dipole matrix element whereas the bare ma- 
trix element is exhibited by a three-pronged vertex. The 
Coulomb interaction matrices (njljEl\\V\\E'l'nili) and 
(E'l' El\\V\\niliTijlj) describe the so-called time-forward 
and time-reverse correlation processes which arc exhib- 
ited by the second and third diagrams (from left to right). 
In the time-forward process, the photon absorption is 
followed by the inter-electron interaction in the form of 
creation of the virtual electron-hole pair in the neigh- 
bouring sub-shell. In the time-reverse process, the virtual 
electron-hole pair is created before the photon absorption 
takes place. Because the time-forward process is real in 
a sense that it conserves the energy of the system, while 
the time-backward process is virtual, the time-forward 
process makes stronger contribution to the photoioniza- 
tion process. However, for the completeness and gauge 
invariance of the theory, both processes should be taken 
into account. 



To include inter-shell correlation effects, we employ the 
RPA model [16]. In this approximation, the reduced 
dipole matrix clement (5) is replaced by its correlated 
counterpart (El\\D\\nili) which accounts for correlation 
between different valence sub-shells. This correlated ma- 
trix element is found as a solution of the system of the 
integral equations: 

(El\\D\\mli) = (El\\r\\niU) + l lim dE ' ( 7 ) 

n'l' 
njlj 

'(E'l'WDWnjl^injljElWVWE'l'mk) 
u - E' + e nj i J + is 

(njljWDWE'l^ipl'ElWVWnjljmli) ' 
lo + E' - e„ 3 / 3 

Here the combined sum/integral sign incorporates both 
the summation over the discrete excited states n'l' with 
the energy e n q> and the integration over the continuum 
dE' from the threshold to infinity. The Coulomb matrix 
contains both the direct and the exchange parts V = 




FIG. 1: Graphical representation of the RPA equations (7). 
Left: non-correlated dipole matrix element. Center: time- 
forward process. Right: time-reverse process. 

We solve the system of integral equations (7) using 
a slightly modified version of the computer code [22]. 
While in the original code correlations between two sub- 
shells could only be taken in account, the present version 
allows to include an arbitrary number of sub-shells. The 
energy integration in the time-forward term of Eq. (7) 
(second line) contains a pole and the RPA matrix ele- 
ment acquires an imaginary part and therefore an extra 
phase arg(£ , Z||_D||njZ i ) . However, this phase does not en- 
ter the partial photoionization cross section n^Z,; — > El 
which is obtained from the squared matrix element. To 
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get access to the phase information, one has to evaluate 
the angular asymmetry parameter /? which contains the 
phase difference between the two photoionization chan- 
nels I = l{± 1 when li ^ [16]. The photoelectron group 
delay, which is the energy derivative of the phase of the 
complex photoionization amplitude, gives an alternative 
access to the phase information. It is evaluated as 



dE 



argf(E) = Im f(E)/f(E) 



(9) 



Here the photoionization amplitude f(E) is given the 
partial wave expansion 



f(E) ex e iS H- l Y lm (k)(-l) m 
x (El\\D\\ml, 
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-m TO; 



(10) 



The amplitude f{E) is evaluated in the forward direc- 
tion k\\z, which is usually the case in the attosecond 
time delay measurements. In this case, Y) m (fc||z) = 
Z(47r) _1 / 2 i5 m o and hence m, = also. It has to be noted 
that the phase of the amplitude (10) contains the con- 
tribution of the HF phases Si in both photoionization 
channels I = li ± 1 as well as the RPA correction due 
to the imaginary part of the RPA dipolc matrix element 
(El\\D\\riili) . Thus the associated group delay is labeled 
HF+RPA when the numerical results are presented in 
the following section. 



III. NUMERICAL RESULTS 
A. Neon 2s and 2p sub-shells 

On the top panel of Fig. 2 we present the partial 
photoionization cross-sections of valence shell photoion- 
ization of Ne. The HF cross-sections are shown by the 
dashed (blue) lines and the RPA cross-sections are ex- 
hibited by the solid (red) line. The recommended ex- 
perimental data by Bizau and Wuillcumier [23] are dis- 
played with error bars. In the RPA calculation, we sub- 
stitute the HF bound state energies with the the ex- 
perimental ionization thresholds £2p 3/ , = 21.56 eV and 
( -2s = 48.47 eV [24] which are indicated on the upper 
boundary of the panel. We see that account for the RPA 
correlation between the 2s and 2p sub-shells improves the 
calculated cross-sections and makes then closer to the ex- 
perimental data. 

We note that even though agreement between the- 
ory and experiment is improved in the RPA model, 
there is a visible difference between the calculated and 
measured cross-sections, especially for the 2s sub-shell. 
This difference may arise from the fact that not all the 
many-electron correlations are accounted for by the RPA 
model which includes pairwise electron-hole virtual exci- 
tations. Other processes like admixture of the two-holc- 
one-electron states to the pure one-hole state in the singly 



charged ion are not included in the RPA model. This ad- 
mixture is responsible for the shift in atomic ionization 
potentials relative to the corresponding HF binding en- 
ergies as well as appearance of the satellite lines in the 
photoionization spectra [25]. These effects cannot be ac- 
counted for ab initio in the RPA model. Phenomcnolog- 
ically, they are partly compensated by using the exper- 
imental ionization potentials instead of the HF energies 
Ei in the RPA equations (7). 
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FIG. 2: (Color online) Top: the partial photoionization cross- 
sections of the 2s and 2p sub-shells of Ne. The HF and RPA 
calculations are shown by the dashed (blue) and solid (red) 
lines, respectively. The recommended experimental data by 
Bizau and Wuilleumier [23] are displayed with error bars. 
Middle: elastic scattering phases in the field of the Ne + ion 
for the 2s — > Ep and the dominant 2p — > Ed channels (dotted 
blue line) and the RPA phases (solid red line). The thin dot- 
ted line visualizes the Coulomb phase with Z = 1. Bottom: 
the phase derivatives are converted to the units of the group 
delay. The vertical bar at the photon energy of 106 eV visu- 
alizes the relative time delay between the 2p and 2s sub-shells 
of 21 as as measured by Schultze et al [1] 

On the middle panel, we show the elastic scattering 
phases in the field of the Ne + ion for the 2s —> Ep and the 
dominant 2p —> Ed channels. The HF phases S p (E) and 
5d{E) are plotted with the dashed (blue) line. The RPA 
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phases arg(fcp||_D||2s) and arg(fcd||£)||2p) are displayed 
with the solid (red) line. The thin dotted line visualizes 

the Coulomb phase <Ji{E) = argT (l + l- iZ eS /V2E\ 

with the effective charge Z e g = 1. This phase shows the 
contribution of the long-range Coulomb potential to the 
HF phase 61(E) which is strongly dominant at small ki- 
netic energies of the photoelectron. The phase shift due 
to the short range potential, i.e. the difference of the total 
phase and the Coulomb phase, is related to the quan- 
tum defect according to the Levinson-Seaton theorem 
Si(k — s- 0) — cri(k —> 0) = /i;(oo)7r [26]. For a neutral tar- 
get, the scattering phase at zero energy is related to the 
number of the bound target states Ni by the Levinson's 
theorem 8i(k — > 0) = Nitt. In the absence of the Coulomb 
potential, the 2s — > Ep phase would tend to one unit of n 
at k — > as there is one occupied np sub-shell in the Ne + 
ion with n = 2. With the Coulomb potential taken into 
account, S2 S ^kp(k —> 0) — (X2 S ^kp(k — > 0) = 0.88-7T where 
/i;=i = 0.88 is the quantum defect calculated from fitting 
the np orbital energies in the Ne + ion e np oc — (n — /i; = i) 2 
for n > 2. As the Coulomb phase tends to zero rapidly 
away from the threshold, the HF phase stays rather flat 
at the value determined by the corresponding quantum 
defect. We may associate this behavior with the Levinson 
theorem, even though this theorem is strictly valid only 
at k — >• 0. Similarly, the the 2p — > Ed phase would tend 
to zero as there are no occupied d-shells left behind. The 
Coulomb logarithmic singularity changes this behaviour 
radically and sends the scattering phases to large nega- 
tive values near the threshold. It has to be noted that 
when the Coulomb behavior of the phases and associated 
group delays becomes dominant at low photoelectron en- 
ergies, the measurement-induced components in the ex- 
periments, i.e. the Coulomb-laser coupling in attosecond 
streaking [12, 13] or the continuum-continuum contribu- 
tion in interfcromctric two-photon measurements [2, 5, 6], 
will be large and the group delay as presented in this pa- 
per can only be accessed if those corrections are properly 
accounted for. 

The RPA phase in the 2p — > Ed channel is hardly 
distinguishable from zero. This observation is consistent 
with a very small change that the RPA correction causes 
to the partial photoionization cross-section shown on the 
top panel. The RPA phase in the 2s — > Ep channel is 
large but rather flat and changes slowly with the photon 
energy. This is consistent with the 2s partial photoioniza- 
tion cross-section which is affected by by the inter-shell 
correlation with 2p across the whole range of the studied 
photon energies. 

The bottom panel of Fig. 2 displays the photoelec- 
tron group delay calculated as the energy derivative of 
the phase of the photoionization matrix element. The 
HF group delay in the dominant photoionization chan- 
nel is calculated as t^ f (as) = k~ 1 d8i/dk x 24.2. Here 
E = k' 2 /2 is the photoelectron energy in atomic units 
and one unit of time is equal approximately to 24.2 as. 
In the existing code, the continuous electron orbitals are 
calculated on the regular momentum grid and numeri- 



cal differentiation over the momentum, rather than en- 
ergy, is easier to implement. The fine grid of 0.05 au 
of photoelectron momentum is sufficient for an accu- 
rate numerical differentiation. Similarly, the combined 
RPA+HF time delay is calculated as T^ PA+HF (as) = 

Im fc _1 /'(fc)//(fc) x 24.2 . Here the photoionization am- 
plitude (10) is evaluated in the z-axis direction. We see 
that the HF time delay in the dominant 2p — > Ed channel 
accounts for almost the whole time delay in photoemis- 
sion from the 2p sub-shell. There is some oscillation visi- 
ble due to the autoionizing resonances near the 2s thresh- 
old which is absent in the HF approximation. Overall, 
the 2p time delay is always positive and rapidly decreas- 
ing function of the photon energy. This is explained by 
the monotonously decreasing HF phase in the d-partial 
wave which is driven by the Coulomb logarithmic singu- 
larity. The situation is different in the 2s — > Ep channel. 
Here the HF phase crosses over from the Coulomb behav- 
ior at low photoelectron kinetic energy to the Levinson 
behavior at larger energies. In result, the phase deriva- 
tive and, consequently, the time delay change their sign 
from positive and negative towards the larger photon en- 
ergies. The RPA correction to the time delay is always 
negative. Hence the photoemission from the 2s sub-shell 
seems to be ahead of that of the 2p sub-shell at around 
100 eV photon energy mark where the measurement of 
Schultze et al [1] was taken (shown as a vertical bar in the 
figure). According to (2), to make a comparison of the 
present calculation with the experiment, we have to add 
to the Wigner time delay difference between the 2p and 
2s sub-shells Ar\y = 8.4 as with the difference between 
the corresponding CLC corrections Atclc = 3.5 as [14]. 
The resulting time delay difference Ar s = 11.9 as which 
is very similar to that reported in [10] by only half of the 
experimental value of 21 ± 5 as. 



1. Argon 3s and 3p sub-shells 

An analogous set of data for Ar 3s and 3p sub-shells 
is shown in Fig. 3. On the top panel we make a com- 
parison of the HF ( dashed blue line) and the RPA (solid 
red line) partial photoionization cross-sections with the 
experimental data by Mdbus et al [27] for 3s sub-shell 
and by Samson and Stolte [28] for the sum of 3s and 
3p sub-shells. The experimental ionization thresholds 
e 3p3/2 = 15.76 eV and e 3s = 29.24 eV [24] are indicated 
on the upper boundary of the panel. These partial pho- 
toionization cross-sections are qualitatively different from 
those of Ne shown in Fig. 2. Firstly, the 3p cross-section 
in Ar displays the Cooper minimum whereas the nodeless 
2p orbital does not [29]. Second, the inter-shell correla- 
tion changes completely the 3s cross-section which also 
display a deep Cooper-like minimum at a slightly smaller 
photon energy. The RPA calculation reproduces these 
features in fair agreement with the experiment. 

The HF phases in Ar behave similarly to the analo- 
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gous case of Ne except that the 3s — s- Ep phase would 
tend to 2-7T in the absence of the Coulomb singularity as 
there are two occupied np-shells in the Ar + ion. With 
the Coulomb potential taken into account, 53 S ->kp(k — > 
0) — (T3s_^fc p (fc — > 0) = 1.737T where the corresponding 
value of the quantum defect in Ar + is /iz=i = 1.73. The 
RPA phases in Ar are very different from Ne. When 
the cross-section goes through the Cooper minimum, the 
corresponding phase makes a jump of tt in the 3s — > Ep 
amplitude, and — tt in the 3p — > Ed amplitude. This 
jump is easy to understand. If the amplitude was real 
and had a node, it would simply change its sign which 
would amount to adding a phase factor of tt in the com- 
plex number representation. Incidentally, this jump was 
investigated in an earlier model calculation [30] which es- 
tablished validity of the attosecond streaking technique 
for the phase measurements. 

This jump of tt has a dramatic effect on the time delay 
which is shown on the bottom panel of Fig. 3. It drives 
the time delay in the 3s sub-shell to very larger numbers 
in several hundreds of attoscconds. The situation is less 
dramatic for the 3p sub-shell. Here the normally weak 
transition Zp —> Es takes over near the Cooper minimum 
of the strong 3p — > Ed transition and the resulting time 
delay docs not go below —100 as. We note that there is 
a strong variation of phase near the autoionization res- 
onances in the 3p photoionization which is seen on the 
top panel of Fig. 3. We do not show this variation in 
the phase and time delay plots for clarity of presenta- 
tion. Anyway, this resonances are far too narrow to be 
detected in time delay measurements at present energy 
resolution. 

One can compare significant time delay near the 
Cooper minimum with the delay time in Breit-Wigner 
resonant scattering td = 2/Y with V being the resonant 
width at half maximum of the cross-section [31]. In the 
case of the Cooper minimum in Ar, which is roughly 
0.5 a.u. of energy wide, the time delay is expected to 
be 4 atomic units of time which equates to about 100 as. 
Of course, this is a very rough estimate and the actual 
time delay is not constant but varies across the Cooper 
minimum. The steepness of this variation can only be 
estimated from an accurate numerical calculation. 

On the upper boundary of the bottom panel, we in- 
dicate the photon energies corresponding to the SB 22 
to 26 of the titanium:sapphire laser at 800 nm used in 
the two-photon interferometric experiments [2, 6] We see 
that at this photon energy range, the RPA correction 
changes completely the sign of the relative 3p/3s time 
delay. In the HF approximation, the 3p photocmission 
is delayed more that the 3s ones. The inter-shell corre- 
lation changes this ordering completely With the RPA 
correction, it is the 3s that is delayed more than the 
3p. This is an important, strong and qualitative result 
which is related to the Cooper minima in the correspond- 
ing partial photoionization cross-sections. This result is 
confirmed by an alternative time-independent calculation 
by Dahlstrom et al. [10] with a similar account for many- 



electron correlations as in RPA. As compared to the orig- 
inal calculation presented in [10], the group delay data 
shown on the bottom panel of Fig. 3 are corrected for 
the experimental ionization potentials [32]. Without this 
correction, the HF ionization potential of the 3s sub-shell 
e-js = 34.7 eV makes the SB 22 inaccessible. 
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FIG. 3: (Color online) Top: the partial photoionization cross- 
sections of the 3s and 3p sub-shells of Ar. The HF and RPA 
calculations are shown by the dashed (blue) and solid (red) 
lines, respectively. The experimental data for 3s [27] and for 
3s + 3p [28] are displayed with error bars. Middle: elastic scat- 
tering phases in the field of the Ar + ion for the 3s —¥ Ep and 
the dominant 3p — > Ed channels (dotted blue line) and the 
RPA phases (solid red line). Bottom: the phase derivatives 
are converted to the units of the group delay. The green as- 
terisks display the calculation [10] corrected for experimental 
ionization thresholds. 



A strong modification of the relative time delay be- 
tween the 3p and 3s sub-shells in Ar is more clearly 
seen in Table I where we present the time delay dif- 
ference T3 S — T3 P in the HF and RPA approximations 
and compare it with the experimental data of Guenot 
et al. [6]. Even a fairly large uncertainty of ±50 as 
cannot reconcile the experimental data with neither of 
the calculations. In the same table, we present results 
of a multi-configurational Hartrcc-Fock (MCHF) close- 
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SB 22 24 26 

to (eV) 34.1 37.2 40.3 



T 3a 



T w ( as ) 








HF 


3 


-36 


-38 


RPA 


76 


53 


215 


MCHF 


45 


10 


-5 


Expt 


70 


-30 


50 



TABLE I: Relative time delay between the photoemission 
from the 3s and 3p sub-shells Atw = Tw — in Ar at 
three fixed photon energies corresponding to the SB 22 to 26 
in the experiment of Guenot et al. [6]. The experimental un- 
certainty is ±50 as. The MCHF calculation [33] has typical 
error bars of ±40 as due to pseudo-resonance structure. 



coupling calculation [33]. In this calculation, the Cooper 
minimum was displaced to significantly larger photon en- 
ergies which were not probed experimentally. Hence the 
time delay difference at the SB 26 was not affected by this 
minimum as strongly as in the present RPA calculation. 



2. Krypton 4p, 4s and 3d sub-shells 

Our results for the 4p, 4s and 3d photoionization of 
Kr are displayed in Fig. 4. On the top panel we 
make a comparison of the HF ( dashed blue line) and 
the RPA (solid red line) partial photoionization cross- 
sections with the experimental data of Ehresmann et al. 
[34] for 4s and of Samson and Stolte [28] for Ap ± 3d 
(error bars). The data from Aksela et al. [35] for 3d 
are displayed with asterisks. The experimental ioniza- 
tion thresholds e 4p3/2 = 14.00 eV, e 4s = 27.51 eV [24] 
and e3d E , 2 = 93.83 eV [36] are indicated on the upper 
boundary of the panel. The 4p and 4s cross-sections in 
Kr behave similarly to the 3p and 3s cross-sections in Ar 
(see the top panel of Fig. 3). The Ap — > Ed cross-section 
goes through its Cooper minimum which is offset some- 
what by the weaker 4p — > Es channel. So the total Ap 
cross-section displays a shoulder rather than a true mini- 
mum. The 4s cross-section is driven strongly by its inter- 
shell correlation with 4p to a very deep minimum which 
is missed completely in the HF approximation. The 3d 
cross-section from its threshold displays a strong maxi- 
mum associated with its shape resonance. This resonance 
is known to be due to electron correlation within a single 
shell [37] and indeed the 3d photoionization cross-section 
is well described by the HF approximation. 

The HF phases in Kr (middle panel of Fig. 4) behave 
similarly to the analogous cases of Ne and Ar except that 
the 4s — > Ep phase would tend to 37r and the Ap — > Ed 
phase would tend to ir in the absence of the Coulomb 
potential. With this potential, the HF phases are de- 
termined by the corresponding quantum defect values 
= 2.67 and ^/ =2 = 1-04. The RPA phases in Kr 
are also similar to Ar. Every time the cross-section goes 
through the Cooper minimum, the corresponding phase 



makes a jump of it: upwards in the 4s — > Ep amplitude 
and downwards in the 4p — > Ed amplitude. The RPA 
phase in the 3d —> Ef transition is rather stationary. 
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FIG. 4: (Color online) Top: the partial photoionization cross- 
sections of Kr. The HF and RPA calculations are shown 
by the dashed (blue) and solid (red for 4s, 4p and green for 
3d)) lines, respectively. The experimental data by Ehresmann 
et al. [34] for 4s and by Samson and Stolte [28] for 4p ± 3d 
are displayed with error bars. The data from Aksela et al. 
[35] for 3d are displayed with asterisks. Middle: elastic scat- 
tering phases in the field of the Ar + ion for the 4s — > Ep and 
the dominant 4p — > Ed channels (dotted blue line) and the 
RPA phases (solid red line). Bottom: the phase derivatives 
are converted to the units of the group delay. 

This behavior of the phases translates into the corre- 
sponding time delays plotted on the bottom panel of Fig. 
4. The RPA time delay in 4p sub-shell is not dramatically 
different from the HF calculation. Even though the dom- 
inant 4p — > Ed transition displays a Cooper minimum, 
it is offset by the weak Ap — > Es transition and is not as 
prominent in the total Ap cross-section as in the 3p cross- 
section of Ar. There are some variation of the time delay 
near the autoionizing resonances close to the 4s thresh- 
old which are seen in the RPA calculation but not in HF 
one. The time delay in the 3d sub-shell is almost entirely 
due to intra-shcll effects and the HF and RPA results 
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are very close. The situation is very different in the 4s 
sub-shell where the time delay is strongly affected by the 
inter-shell correlation with the Ap sub-shell and reaches 
300 as in its peak. Similarly to Ar, there is a complete 
reversal of the relative time delay between the Ap and 4s 
sub-shells in the RPA calculation in comparison with the 
HF one. 



3. Xenon 5p, 5s and Ad sub-shells 

The analogous set of data for the bp, 5s and Ad sub- 
shells of Xe is presented in Fig. 5. On the top panel we 
compare the partial photoionization cross-sections in the 
HF (dashed blue line) and RPA ( solid red line) approx- 
imations with the experimental data [38, 39] which are 
shown with the blue asterisks for 5 s and error bars for bp 
and Ad. The experimental ionization thresholds es P3/2 = 
12.13 cV , e 2s = 23.40 eV [24] and e 4 d 5/2 = 67.50 eV [40] 
are indicated on the upper boundary of the panel. 

Below the Ad ionization threshold, the 5s and bp cross- 
sections in Xe behave similarly to to the 4s and Ap sub- 
shells in Kr (top panel of Fig. 4). However, above this 
threshold, the Ad sub-shell goes through a very steep 
shape resonance, sometimes even called a "giant reso- 
nance" . This resonance is then turns into a Cooper 
minimum. By strong inter-shell interaction, this behav- 
ior is replicated in the bp and 5s partial photoioniza- 
tion cross-sections which are well reproduced by the RPA 
calculation. Accordingly, the corresponding RPA phases 
displays steep tt jumps (middle panel) which are reflected 
in the corresponding time delays (bottom panel) . In the 
case of the 5s sub-shell, the RPA phase jump near the 
Cooper minimum mergers with the Coulomb singularity 
and produces a very large, nearly 300 as time delay at the 
photon energies below 30 eV. The bp sub-shell shows a 
large and negative time delay due to its Cooper minimum 
at around 50 eV. Both the 5s and hp sub-shells display a 
large and negative time delay near the local cross-section 
minima around 150 eV induced by the correlation with 
the Ad sub-shell. The time delay in the Ad sub-shell is 
driven from the strongly positive due to the Coulomb sin- 
gularity at low photon energies to a large negative jump 
near the Cooper minimum at about 180 eV. At larger 
energies, the cross-sections are rather structureless and 
there is no significant time delay variations. 

A phase jump of 7r, smoothed by the interaction 
between the two channels, has already been observed 
both theoretically and experimentally by analyzing the 
anisotropy parameter in photoionization of Xe bp sub- 
shell [41]. This parameter contains the phase shift be- 
tween the two photoionization channels with I = ± 1. 
In the case of bp photoionization, these are bp — > Ed 
and bp — > Es transitions. Their partial photoionization 
cross-sections and the relative phase shift are presented 
on the top and bottom panels of Fig. 6. On both panels, 
we show the present RPA and HF calculations displayed 
with the solid red and blue dotted lines, respectively. On 
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FIG. 5: (Color online) Top: the partial photoionization cross- 
sections of Xe. The HF and RPA calculations are shown by 
the dashed (blue) and solid (red for 5s, 5p and green for Ad) 
lines, respectively. The experimental data from Becker et al. 
[38] and Fahlman et al. [39] are shown with asterisks for 5s 
and error bars for 5p and Ad. Middle: elastic scattering phases 
in the field of the Ar + ion for the 5s — > Ep and the dominant 
5p —¥ Ed and Ad — > Ef channels (dotted blue line) and the 
RPA phases (solid red line). Bottom: the phase derivatives 
are converted to the units of the group delay. 



the bottom panel, we exhibit the RPA (open circles) and 
HF (filled circles) phase shifts reported by [41]. 

On the top panel of Fig. 6 we observe a significant 
shift of the Cooper minimum in the bp — > Ed channel 
towards the lower photon energies and appearance of the 
secondary minimum due to the correlation with the Ad 
sub-shell. In the meantime, the inter-shell correlation 
does not change the bp — > Es partial photoionization 
cross-section in such a dramatic way. Accordingly, on 
the bottom panel of Fig. 6, we see a strong variation of 
the RPA phase shift with the two successive jumps near 
the Cooper minima of the bp — > Ed cross-section. In the 
meantime, the HF calculation returns quite a smooth and 
monotonous phase shift. Agreement between the two sets 
of calculations, the present and the one reported by [41], 
is rather good. A small shift between the present calcula- 
tion and the reference one is most likely due to scanning 
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FIG. 6: (Color online) Top: Partial photoionization cross- 
sections of Xe in the 5p —¥ Ed and 5p — > Es channels in the 
RPA (solid red line) and HF (dotted blue line) approxima- 
tions. Bottom: Phase shift between the partial 5p —5- Ed and 
5p — > Es waves. The present RPA and HF calculations (solid 
red and blue dotted lines, respectively) are compared with the 
RPA and HF calculations reported in Zimmermann et al. [41] 
(open and filled circles, respectively). 



and digitizing the analog data of Fig. 3 in Ref. [41] . 



IV. CONCLUSION 

In the present work, we perform a systematic study of 
the photoemission time delay from the valence shells of 
noble gas atoms in sequence from Ne to Xe. We cover 
the photon energy range from the ionization threshold 
to 200 cV. We test the accuracy of our calculation by 
making comparison with available partial photoioniza- 
tion cross-sections. We derive the complex phase of the 
photoionization amplitude in the non-relativistic HF and 
RPA calculations and convert it to the photoelectron 
group delay by taking the energy derivative. 

The time delay results display a very diverse land- 
scape due to an interplay of three major factors. The 
first two are the logarithmic Coulomb singularity and the 
Lcvinson theorem which drive the photoelectron scatter- 
ing phase in the field of the singly charged ion. The 
third factor is the phase jump of it near the Cooper min- 
imum which is smoothed by the inter-shell interaction. 
The two former factors are revealed in the HF calcula- 
tions whereas the third one is most vividly reflected in 
the RPA calculations. Experimentally, photoionization 
measurements near the Cooper minima may be challeng- 
ing but it is the area where the time delay effects are 
expected to be largest. 



These time delay results are compared with experi- 
mental data derived from the attosecond streaking mea- 
surement [1] and the two-photon interferometric tech- 
nique [6]. This comparison is inconclusive as the dif- 
ference between the theoretical and experimental results 
clearly exceeds the reported error bars. We are fairly 
confident about the accuracy of the present calculation 
which is tested by comparison of the partial photoion- 
ization cross-sections with a large set of independent ex- 
perimental data and the angular asymmetry parameters 
as in the case of Xe [41]. It is hard to give a numerical 
estimate on the accuracy of the group delay results. In 
lighter atoms we expect it to be within 10%. For Xe, it 
may be more significant as suggested by larger difference 
between the calculated and experimental cross-sections. 
Even for this heaviest of the atoms studied in the present 
work, the relativistic effects are not expected to change 
considerably the complex phase [41] and hence the asso- 
ciated group delay. It is therefore an open question why 
the time delay results cannot be verified experimentally 
even after the corresponding CLC or CC corrections are 
made. Such a verification would be a very welcoming de- 
velopment both for the attosecond time delay measuring 
techniques and the complete theory of atomic photoion- 
ization. 
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